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Abstract 

A transfer matrix function representation of the fundamental solu- 
tion of the general-type discrete Dirac system, corresponding to rect- 
angular Schur coefficients and Weyl functions, is obtained. Connec- 
tions with Szego recurrence, Schur coefficients and structured matrices 
are treated. Borg-Marchenko-type uniqueness theorem is derived. In- 
verse problems on the interval and semiaxis are solved. 
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1 Introduction 

In this paper we deal with a discrete Dirac-type (or simply Dirac) system: 
y k +i(z) = (I m + izjC k )y k (z) (k e N ) , (1.1) 

where No stands for the set of non-negative integer numbers, I m is the mxm 
identity matrix, "i" is the imaginary unit (i 2 = — 1) and the mxm matrices 
{C k } are positive and j-unitary: 



C k > 0, CkjCk = j, j 



I mi 

~Im 2 J 



(mi + m 2 = m, m 1; m 2 ^ 0). 

(1.2) 
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Discrete systems are of great interest and their study is sometimes more 
complicated than the study of the corresponding continuous systems (see, 
e.g., [TH3|[7t[T2] and references therein). The subcase mi = m 2 of system (11. ip 
(satisfying (II .2p ) corresponds to the self-adjoint Dirac-type systems, which 
were studied in [14] (and the subcase j = I m of system ( 11. ip corresponds to 
the skew- self-adjoint Dirac-type systems, an important subclass of which was 
investigated in The analogies between system (II. ip and continuous 

Dirac-type systems are also discussed in [TU[20j|23] m detail. Here we follow 
the paper [15] on the continuous case, where mi does not necessarily equals 
m 2 and the m 2 x mi Weyl matrix functions are, correspondingly, rectangular. 

It is essential that Dirac system ( II. ip . ( II ,2p is equivalent to the very 
well-known Szego recurrence (see, e.g., [101128]). This connection is discussed 
in detail in Section [2j Inverse problems for the subcase of the scalar Schur 
(or Verblunsky) coefficients were studied, for instance, in 0128] (see also 
various references therein), and here we deal with the rectangular matrix 
Schur coefficients. 

In this paper Im denotes image of a matrix (or operator), a (A) stands 
for the spectrum of A and "span" stands for the linear span. 

2 Dirac system and Szego recurrence 

The next simple proposition is essential for our future research and could be 
of independent interest in the theory of functions (and powers, in particular) 
of matrices, which is developed in a series of works (see, e.g., [6l|29] and 
references therein). 

Proposition 2.1 Let an m x m matrix C satisfy relations 



0, CjC ./ (./ ./• j '). 



(2.1) 



Then the following relations hold for all sGl: 



C s > 0, C s jC s = j. 



(2.2) 



Proof. Since C > 0, it admits a representation 



C = u*Du, 



(2.3) 
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where D is a diagonal matrix and 

D > 0, u*u = uu* = I m . (2.4) 
We substitute f 12 . 3 It into the second equality in ( 12.1 ft to derive 

u*Duju*Du = j, 

or, equivalent ly, 

DJD = J, J=J* = J' 1 := uju*. (2.5) 

Formula (I2.5P yields D^ 1 = JDJ and, taking power s of the both parts of 
this equality, we obtain 

D~ s = JD S J, D S JD S = J. (2.6) 

Finally, using (12.411 — (12.6ji we have 

u*D s uju*D s u = j. (2.7) 

■ 

We substitute s — 1/2 and apply Proposition 12.11 to matrices Ck in order to 
obtain the next proposition. 

Proposition 2.2 Let matrices Ck satisfy ( 11. 21) . Then they admit represen- 
tations 

C h = 2P(k)*P(k)-j, 0(k)jP(k)* = I mi , (2-8) 
C fc =j + 2 7 (fcr 7 (fc), 7(AOi70r = -U, (2-9) 

where (3(k) and j(k) are mi x m and x m matrices given by (12. lQf) and 
(12. lip , respectively. 

Proof. We note that matrices Ck satisfy conditions of Proposition 12.11 and 
so (TjLD holds for C = C k . Next we put 

P(k):=[l mi 0}Ct /2 (2.10) 
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yl/2 
J k ■ 



and take into account the equality 

C k = Cl /2 (2[l mi 0]*[/ mi 0}-j)cl 

Now, representation (12. 8p is apparent from (I2.2p taken with s — 1/2. In a 
similar way, formula ( 12. 2 p and equality I m = j + 2 [0 J TO2 ] [0 J m2 ] imply 
representation (12.91) for 

7(*)=[0 Im 2 ]C l k 12 . (2.11) 



Now, we will consider interrelations between Dirac system ( II. ip . ( 11.21) 
and Szego recurrence, which is given by the formula 



X k+ \(\) — T> k H k 



A/ mi 







L m,2 



X k {\), 



(2.12) 



where 



Hi, 



'mi Pk 

* T 
Pk 1 rn2 



V k = diag|(/ mi -p fc p^) 5 , (l m2 -p* k p k ) 2 |, (2.13) 



and the mx x m 2 matrices p k are strictly contractive, that is, 

llPfcll < !■ 



(2.14) 



Remark 2.3 When m x = m 2 = 1, one easily removes the factor (1 — 
|pfc| 2 ) _1//2 2n ( I2.12p to obtain systems as in J^lSi, where direct and inverse 
problems for the case of scalar strictly pseudo-exponential potentials have 
been treated. The square matrix version (i.e., the version where mi = m 2 ) 
of Szego recurrence, its connections with Schur coefficients and applications 
are discussed in /HG^ (see also references therein). For the rectangular ma- 
trices p k see, for instance, [Wj . We note that V k H k is the so called Halmos 
extension of p k (see fTffi p. 167J), and that the matrices D k and H k commute 
(which easily follows, e.g., from JTU, Lemma 1.1.12]). The matrix V k H k is 
j -unitary and positive, that is, 

V k H kJ H k V k = H k V kJ V k H k = j, (2.15) 
V k H k > 0. (2.16) 



4 



According to [TTJ Theorem 1.2], any j-unitary matrix C admits a representa- 
tion, which is close to Halmos extension. More precisely, partitioning C into 
blocks C = {cik}f k= i we see that the mi x mi block en and the m 2 x 777.2 
block C22 are invertible. Then, putting 

P = C12C22 1 = (cn)*C2i, ui = (J mi - pp*) 1/2 cn, u 2 = (/ m2 - p*p) 1/2 c 22 , 
we have the respresentation: 



C = VH 



ui ' 

u 2 



L m,i 1 



H 



X? = diag{(/ mi -pp*) \ (J m2 -p*p) h ], 



imi P 

P I'nri2 

p*p < J m2 . 



(2.17) 
(2.18) 



Although relations ( I2.15p - (l2.17j) are well-known, we could not find in the 
literature a statement, which is converse to (I2.15p . ( I2.16p . Hence, we prove 
it below. 

Proposition 2.4 Let an m x m matrix C be j-unitary and positive. Then 
it admits a representation 



C = VH, 



(2.19) 



where H and V are of the form (12.171) and (I2.18P {i.e., the last factor on the 
right-hand side of the first equality in (I2.17P is removed) . 

Proof. Recall that C admits representation (12.171) . We fix a unitary matrix 
U such that T>H = UDU*, where D is a diagonal matrix, D > 0. Then, 
relations C = C* and (12.171) yield the equality 



Ml 


0" 




u\ 


0" 





u 2 _ 







U* 2 _ 



UDU* 

which we rewrite in the form 

DU = U*D, U := U* 



UDU* 



Ui 
u 2 



u . 



(2.20) 
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According to (I2.20p . DU is a selfadjoint matrix, and so D 1 I 2 UD 1//2 is a 
selfadjoint matrix too, that is, there is a representation 



D 1/2 UD- 1/2 = UDJJ*, 



(2.21) 



where U and D\ = D* are unitary and diagonal matrices, respectively. The 
definition of U in ( I2.20p implies that U is unitary. Therefore, in view of 
f !2.2ip . Di is linear similar to a unitary matrix, that is, its entries are ±1. 
Moreover D x > 0, since C > and formulas (BTTJ) . ff^"20|) and (1231]) yield 



«i ' 
u 2 



UDUU* = UD l/2 tjD 1 U*D ll2 U* 



(2.22) 



From the inequality D\ > and the fact that the entries of equal either 1 
or —1, we have Di = I m . Thus, the last equality in (I2.22p implies C = U DU*, 
that is, f l2~T9|) holds. 



Proposition 12.41 completes Propositions 12.11 and 12.21 on representations and 
properties of Ck- Taking into account (I2.15p . (I2.16P and Proposition \2A\ we 
rewrite Szego recurrence f)2.12p in an equivalent form 



Xk+i(\) — Ck 



A/ mi 







'm 2 



X k (X), k e N c 



C k > 0, C k jC k = j. 



(2.23) 
(2.24) 



Using (I2.24p we see that the matrix functions Uk, which are given by the 
equalities 



U :=I m , U k+1 :=iU k C k j = Y[(iC r j) (k>0) 

r=0 

are also j-unitary. From (12.241) and (I2.25j) we have 

(i + z)U k+1 (I m + izj)C k 
= I m + izU k+1 jU^ r 



(2.25) 



— / n 

I m2 



(I m + izj)- 1 U k 1 



(2.26) 



6 



In view of (12.261) . the function y k of the form 



y k (z) = (i + z) k U k (I m + izj)X k [j^fj (2-27) 

satisfies (II. ip . where yo(z) = (I m + izj)X (z) and C k = jU k+ ijU^} x . Since 
U k+ i is j-unitary, we rewrite C k as 

C k = jU k+1 U* k+lJ , (2.28) 

and so (II. 2p holds. Because of ( I2.25p . (I2.28P and j-unitarity of U k , we have 
jU* k C k U k j = Cl that is, 

C k = (jU* k C k U k j) 1/2 . (2.29) 
The following theorem describes interconnections between systems (II. ip and 

Theorem 2.5 Dirac systems fll.ip . f 1 1 . 2 p and Szego recurrences (I2.23p . (I2.24p 

are equivalent. The transformation DJl : {C k } — > {C k } of Szego recurrence 
into Dirac system, and the transformation of their solutions, are given, re- 
spectively, by formulas (I2.28P and (I2.27p . where matrices {U k } are defined in 
( I2.25p . The mapping OJt is bijective, and the inverse mapping is obtained by 
applying (I2.29P (and substitution of the result into ( I2.25P ) for the successive 
values of k. 

Proof. It is proved already above that the formulas (I2.28P and (I2.27P de- 
scribe a mapping of Szego recurrence and its solution into Dirac system and 
its solution, respectively. Moreover, the mapping 9Jt is injective, since we 
can successively and uniquely recover C k and U k+ i from C k and U k using 
formulas ( I2.29P and (12. 25 p . respectively. 

Next, we prove that 9Jt is surjective. Indeed, given an arbitrary sequence 
{C k } satisfying (ll.2p . let us apply to the matrices from this sequence relation 
( I2.29P (and substitute the result into (I2.25P ) for the successive values of k. 
In this way we construct a sequence {C k }. Since the matrices jUlC k U k j 
are positive and j-unitary, we see, from ( I2.29p and Proposition 12.11 that the 
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matrices C k are also positive and j-unitary. Next, we apply to {C k } the 
mapping 971. Taking into account (I2.25P and (12.291) . we derive 



ju k+1 u; +lJ = ju k c 2 k u* kJ = ju k ( 3 u* k c k u k j)u;j = c k , (2.30) 

that is, DJl maps the constructed sequence {C k } into the initial sequence 
{C k }. Recall that we started from an arbitrary {C k } satisfying ( II. 2p . Hence, 
9Jt is surjective. ■ 



3 Weyl theory: direct problems 

In this section we introduce Weyl functions for matricial discrete Dirac sys- 
tems (II. ip . Next we prove the Weyl function's existence and, moreover, give 
a procedure to construct it (direct problems). Finally, we construct the S- 
node, which corresponds to system ( II. ip . and the transfer matrix function 
representation of the fundamental solution W k . (See, e.g., [2"5T - |2"T] on the 
S- nodes and the transfer matrix functions in Lev Sakhnovich sense.) 

The fundamental m x m solution {TVfc} of (II. ip we normalize by the 
condition 

W (z) = I m . (3.1) 

Similar to the continuous analog of (II. ip in [T5|IT6] (see also canonical system 
case [27J P- 7]), the Weyl functions of system (II. ip on the interval [0, r] (i.e., 
system ( II. ip considered for < k < r) are defined by the Mobius (linear- 
fractional) transformation: 

cp r (z,V)=[0 J m2 ] W r+x {z)- x V{z)([l mi 0}W r+1 (z)- 1 V(z)y\ (3.2) 

where V(z) are nonsingular m x mi matrix functions with property-j. That 
is, V(z) are meromorphic in C + matrix functions such that 

V{z)*V{z) > 0, V(z)*jV(z) > (3.3) 

for all points in C + (excluding, possibly, a discrete set). The first inequal- 
ity in (I3.3p means non-singularity (non- degeneracy) of V and the second 
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inequality is called property- j. Since V is meromorphic, property-j almost 
everywhere in C + and the first inequality in ( 13. 3 p at some z e C + suffice for 
the conditions on V to hold. 

It is apparent from (jl.ip and (13.11) that 



W r+ i(z) = H(I m + izjC k ). (3.4) 

k=0 

In view of (12. 9p and (13.41) we obtain 

W P+1 (i) = (-2r +i n(jW7(*0)- (3-5) 



k=0 



Hence, det W r +i(i) = 0, and we don't consider z = i in this section. 

Remark 3.1 We note that the behavior of Weyl functions in the neighbor- 
hood of z = i is essential for the inverse problems that are dealt with in 
the next section. Therefore, unlike the Weyl disc case (see Notation 3.4), 
in the definition (13.21) of the Weyl functions on the interval we assume that 
V is not only nonsingular with property-j but has also an additional prop- 
erty. Namely, it is well-defined and nonsingular at z = i. We don't use this 
additional property in this section, though, in important cases, it could be 
obtained via multiplication by a scalar function. 

The lemma below shows that transformations tp r (z,V) are well-defined. 

Lemma 3.2 Fix any z e C + such that the inequalities detW r (z) ^ and 
(13.31) hold. Then we have the inequality 

det([/ mi 0] W r+1 (z)- l V(z)^ ^0. (3.6) 

Proof. Using and (J23D we obtain 

(I m + izjC k )*j{I m + izjC k ) = (1 + i(z -z) + \z\ 2 )j + 2i(z - z)j{k)*j{k) 
< (1 -2S(z) + \z\ 2 )j, (1 -2$s(z) + \z\ 2 ) > for z^i. (3.7) 

Since the equality (13 ,4p holds, formula (13. 7p implies that 

(Wr+i{.z)~ 1 )* iW r+ i(z)~ l > (1 - 2%z) + \z\ 2 )- r - l j (z e C+, z + i). 

(3.8) 
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Because of (13. 3p and (13. 8p . we see that V : = W r +\(z)~ x V(z) satisfies the 
inequality V*jV > 0. It is apparent that the same inequality holds for the 
matrix [l mi 0] . In other words, Im W r+ i(z)^ 1 V(z) and Im \l mi 0] are 
maximal j-nonnegative subspaces. Therefore, the inequality (I3.6P follows in 
a standard way from j-theoretic considerations (see, e.g., the proof of (13.45}) 
or the proof of [TU inequality (5.6)] for such considerations). ■ 

Corollary 3.3 The following inequalities hold for the fundamental solution 
W r+1 of flTT]) (where {C k } satisfy flT2])) : 

det W r+1 (z) jt 0, W T+1 {z)- 1 = (1 + z 2 )- r ~ l jW r+1 {zyj (z ± ±i). 

(3.9) 

Proof. Relations (13. 7p and (13 .4p imply that 

W r+1 (z)*jW r+1 (z) = (1 + z 2 ) r+1 j, z = z. 
Hence, using analyticity considerations, we obtain 

W r+l (z)*jW r+1 (z) = (1 + z 2 ) r+1 j, (3.10) 
and (13. 9p is apparent. ■ 

Notation 3.4 The set of values of matrices ip r (z,V), which are given by the 
transformation (13. 2 p where parameter matrices V(z) satisfy f 13 . 3 1) . is denoted 
byj\f(r,z) (or, sometimes, simply M(r)). 

Usually, Af(r, z) is called the Weyl disk. 

Corollary 3.5 The sets J\f(r,z) are embedded (i.e., J\f(r,z) C J\f(r — l,z)) 
for all r > and z G C +; z ^ i. Moreover, for all fk (k > 0) we have 

(Pk(z)*ipk(z) < I mi . (3.11) 

Proof. It follows from Corollary 13.31 that the matrices W r+ i(z), W r (z) 
and (J m + izjC r ) are invertible. Hence formulas (13. 3p and (13. 7p imply that 
V := (I TO + izjC r )^ 1 V(z) satisfies (13. 3p . Therefore, we rewrite (I3.2p in the 
form 

<p r (z,V)=[0 I m2 }W r (z)- 1 V(z)^[l mi 0}W r (z)- 1 V(z)y\ (3.12) 
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and see that ip r (z) G Af(r — l,z) (r > 0). Inequality (13. 1 1 j) is obtained for 
the matrices from J\f(0, z) via substitution of r = into (13.121) . ■ 

Weyl functions of system (II. ip on the semiaxis No of non-negative integers 
are defined in a different and more traditional way (in terms of summability) , 
see definition below. We will show also that the definitions of Weyl functions 
on the interval and semiaxis are interrelated. 

Definition 3.6 The Weyl-Titchmarsh (or simply Weyl) function of Dirac 
system (II. ip (which is given on the semiaxis < k < oo and satisfies (II. 2p ) is 
an 777.2 x m i matrix function (p(z) (z G C + ) ; such that the following inequality 
holds: 

oo 

Y,q{z) k [l mi f(z)*] W k (z)*C k W k (z) 

k=0 

q(z) := (I + \z\ 2 )- 1 . (3.14) 
Lemma 3.7 If (p r (z) G Af(r,z), we have the inequality 

r 

J]g(z) fc [/ mi ip r (zY]W k (zYC k W k (z) 

k=0 

X (I m - ip r (z)*(f r (z)). 

Proof. Because of (II. ip and ( II. 2p we have 

W k+1 (z)*jW k+1 (z) = W k (zy(l m - izC k j)j(l m + izjC k )w k (z) 

= q(z)- 1 w k ( z yjw k (z) + \(z - z)w k ( z yc k w k (z). 

(3.16) 

Using (13. 1 p and (13.161) . we derive a summation formula, which is similar to 
the formula for the case that m\ = rri2, see [T4l formula (4.2)]: 

J2q( z ) k w k (zyc k w k (z) = l±^t{j- q ( z y^w r+x (zyjw r+1 (z)). (3.17) 

On the other hand, it follows from ( 13 .2p that 

= W r+l (z)- 1 V(z)(^[l mi 0}W r+1 (z)- 1 V(z)y\ (3.18) 



J rai 

<f(z) 



< OO, 



(3.13) 



J ITll 

cpJz) 



< 



1 + 1*1 



nz 



(3.15) 
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and so formula (13.31) yields 



[I mi <p r (z)*] W r+1 (z)*jW r+1 ( 



ip r (z) 



> 0. 



(3.19) 



Formulas (l3~TTj) and (l3~T9|) imply (13~T5|) . ■ 

Now, we are ready to prove the main direct theorem. 

Theorem 3.8 There is a unique Weyl function of the discrete Dirac system 
(11.11) . which is given on the semi-axis < k < oo and satisfies (11.21) . This 
Weyl function is analytic and non-expansive {i.e., ip*ip < I mi ) in C+. 

Proof. The proof consists of 3 steps. First, we show that there is an analytic 
and non-expansive function 



(3.20) 



r>0 



Next, we show that (poo(z) is a Weyl function. Finally, we prove the unique- 
ness. 

Step 1. This step is similar to the corresponding part of the proof 
of jini Proposition 2.2]. Indeed, from Corollary 13.51 we see that the set of 
functions <p r (z, V) of the form ( 13. 21) is uniformly bounded in C+. So, Montel's 
theorem is applicable and there is an analytic matrix function, which we 
denote by tpoo(z) and which is a uniform limit of some sequence 



(Poo(z) = lim (p r .(z,Vi) (i G N, r { f, 



lim Ti 

i— >oo 



OO 



(3.21) 



on all the bounded and closed subsets of C+. Clearly, ifoo is non-expansive. 
Since ^, the sets Af(r, z) are embedded and equality (I3.18P is valid, it 
follows that the matrix functions 



Vij(z) 



1 m\ 

(p r (Z, V 



(j > 



satisfy relations (13. 3p . Therefore, using (13.211) we derive that (13. 3p holds for 



Vi >0 o(z) := W u+1 (z) 



<Poo{z). 



(3.22) 



12 



which implies that we can substitute V = Pj j00 and r = into (13. 2p to obtain 

(PooWeAffaz). (3.23) 

Since (I3T23D holds for all i G N, we see that fl3T20|) is fulfilled. 

Step 2. Because of (13.201) . the function ip^ satisfies condition of Lemma 
13. 71 Hence, ( 13. 1 5|) holds for any r > and y? r = <^oo, which implies (13. 13j) . 
Therefore, ip^ is a Weyl function. 

Step 3. It is apparent from (I2.8P that 

W k (zyC k W k (z) > W k (z)*(-j)W k (z). (3.24) 
Using (I3.16P we derive also 

q(z) k W k (z)*(-j)W k (z) > qizy^W^izfi-fiW^iz). (3.25) 
Formulas (13. ip . (I3.24p and (I3.25|) yield the basic for Step 3 inequality 

q(z) k W k (zyC k W k (z) > -j. (3.26) 
Therefore, the following equality is immediate for any g G C™ 2 : 





^£*[0 I m Mz) h W k (z)*C k W k (z) 



k=0 



g = oo. 



(3.27) 



It was shown in Step 2 that <p = p^ satisfies (13. 13j) . According to (I3.13P and 
(I3.27p . the dimension of the subspace L G C m of vectors h such that 



J2h*q(z) k W k (z)*C k W k (z)h < oo 



(3.28) 



k=0 



equals mi. Now, suppose that there is a Weyl function cp ^ ip^. Then we 
have 



Im 



<Poo{z) 



C L, Im 



C L. 



Therefore, dimL > mi (for those z, where <p{z) ^ (p^z)) and we arrive at 
a contradiction. ■ 

Finally, let us construct representations of W r+ i (r > 0) via S'-nodes. 
First, recall that matrices {C k } generate via formula (12.1 ip a set {'y(k)} of 
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the m 2 xm matrices 'y(k). Using {^(k)}, we introduce m 2 (r + l) xm matrices 
T r and m 2 (r + 1) x m 2 (r + 1) matrices K r (0 < r < oo): 





7(0) 




x r (0) 


r r := 


7(1) 


; K r := 


x r (l) 




7(r) 




x r (r) 



(3.29) 



Xr{k) := vy(k)j [ 7 (0)* ... 7 (A; - 1)* 7(*0*/2 ... 0] . (3.30) 
It is apparent from (13.291) and (13.301) that the identity 

K r - k; = ir r jr; (3.31) 

holds. The m 2 (r + l) x m 2 (r + l) matrices A r are introduced by the equalities: 

{0 for n > 0, 

(i/2)J m2 for n = 0, (3.32) 
iI mo for n < 0. 



Proposition 3.9 Matrices K r and A r are linear similar. 



K r = E r A r E;\ 



(3.33) 



Moreover, the similarity transformations E r can be constructed so that 



£ r _i 
X r e~ 



(r > 0), E^T^ = $ r , 2 , $ r , 2 



L m 2 



L rri2 J 



(3.34) 



E = e = 7 2 (0), 



(3.35) 



whereTr tP are m 2 (r + l) xm p blocks of T r = \T r> i r rj2 ] andj p (k) are m 2 
blocks ofj(k) = [71 (fe) 72W] ■ 

Proof. It follows from (l3T29|) . ( 13~30|) and ( 13~32|) that 



xm r 



K = A = -(i/2)/ m2 , det 72 (0)^0, 

x r (r) = i [ 7 (r)i7(0)* ... 7 (r)j7(r - 1)* -/J2]. 



(3.36) 
(3.37) 



14 



We see that (13.35|) and (I3.36P imply f !3.33|) for r = 0. Next, we prove (13. 'd'c 
by induction. Assume that K r _i = E r ^iA r ^iE~\ and let E r have the form 
(I3.34p . where dete~ ^ 0. Then we obtain 



EZ l 



(3.38) 



and, in view of (I3.29p . (I3.32p . (I3.34p . (I3.37p . it is necessary and sufficient (for 
( 13~33|) to hold) that 



Im 2 0] 



L rrri2 



(|X4-i ~(i/2)e-] - ie r [l m2 . 
= i 7 (r)j [ 7 (0)* ... -y(r -!)*]. 
We can rewrite (I3.39P in the form 

X r (A._! + (i/2)/ rm2 ) =i 7 (r)j [ 7 (0)* ... 7 (r - 1)*] 



^V-l 

(3.39) 



r-1 



(3.40) 



We partition X r (r > 1) into two m2 x mi and m.2 x (r — l)m2, respectively, 
blocks 



(3.41) 



and we will need also partitions of the matrices A r _i + (i/2)J rm2 and i? r _i, 
which follow (for r > 1) from fl3T32|) and fETMj) : 



(A r _! + (i/2)/ rm2 ) 



0' 

(A r _ 2 - (i/2)/ (r _ 1)m2 ) 





" " 




' 














_ e r-l_ 



(3.42) 



Using (13.41 p and (I3.42[) we see that (I3.40p is equivalent to the relations 

e~ = — 7 (r)j 7 (r — \)*e~_ x for r > 1; (3.43) 
X r =i( T (r)j [ 7 (0)* ... 7(r-ir]Er-i + e r -[/ m2 ... I m2 ] 



x 



(A--2 — (i/2)/( r _i) m2 ) 




-1 



for r > 1. 



(3.44) 
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Hence, if e~ and X r satisfy (I3.43P and (I3.44p . respectively, and dete~ ^ 0, 
the similarity relation (I3.33P holds. The inequalities det e~ ^ are apparent 
(by induction) from (I3.35p . (I3.43P and the inequalities 

det(7(r)j 7 (r-l)*)^0, (3.45) 

and it remains to prove (I3.45p . Indeed, let j(r)jj(r — l)*g = 0, g ^ 0. Then, 
the subspaces Im j(r)* and span j(r — l)*g are j-orthogonal. The second 
equality in f )2.9p (taken for k = r and k — r — 1) implies that these subspaces 
are also j-negative, have zero intersection and have dimensions m,2 and 1, 
respectively. Thus, span ( 7 (r — l)*g U Im j(r)*) is an rri2 + 1-dimensional 
j-negative subspace, which does not exist. Therefore, the relation (I3.45p . and 
so also equality (I3.33p . is proved. 

Formula (13.351) shows that the second equality in (13.341) holds for r = 0. 
Now, we choose X r (for r = 1) and x~ (for r > 1) so that the second equality 
in (13.341) holds in the case that r > 0. Taking into account (I3.38j) . (I3.4ip and 
using induction, we see that this equality is valid when 

X 1 = 7 2 (l)-eT, x- = 72 (r) - e~ - X r $ r _ 2i2 (r > 1). (3.46) 

■ 

We note that inequalities, which are similar to (13. 6p and (I3.45[) . are often 
required in the study of completion problems and Weyl theory. Therefore, 
the next proposition, which is easily proved using the same considerations as 
in the proof of (13.451) . could be of more general interest. 

Proposition 3.10 Let the m x m matrix J satisfy equalities J = J* = J^ 1 
and have mi > positive eigenvalues. Let m x mi matrices d and d satisfy 
inequalities 

■&*$>0, tf*Jtf>0, tfM>0, d*M>0. (3.47) 

Then we have 

detiTJtf^O. (3.48) 
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Let us substitute (I3.33|) into (13.31j) to derive 



E r A r E~ x - (E*) 1 A^E^ = iT rJ r;. (3.49) 

Multiplying both sides of (I3.49P by E~ x and (E*) 1 from the left and right, 
respectively, we obtain the operator identity 

A r S r - S r A* r = m r3 U* r = i($ r ,l®* r ,l - ®r,2$* r , 2 ), (3-50) 

where 

S r := E~ x (E* r y\ n r := E^Tr = [$ r)1 $ r J . (3.51) 



Definition 3.11 The triple of matrices {A r , S r , H r } forms a symmetric S- 
node if the operator (matrix) identity ( I3.50|) holds, S r = S* and det S r ^ 0. 

The transfer matrix function {in Lev Sakhnovich form), which corre- 
sponds to the S-node, is given by the formula 



w A {r, 



\) = i m -iju;s; 1 (A r -xi ir+1)m2 ) %. (3.52) 



Remark 3.12 A symmetric S-node corresponding to Dirac system (11. ip 
(which satisfies ( 11.21) ) on the interval < k < r is constructed using for- 
mulas (I3.32p and (I3.5ip . where T r is given in ( I3.29p . 

Recall that S"-nodes, transfer matrix functions wa and the method of operator 
identities are introduced and studied in [2"4"H2"7] (see also references therein). 
For r > introduce projectors: 

Pi ■= [I rm2 0] , P 2 = P:=[0 ... I m2 ] . (3.53) 

Since E~ x is a block lower triangular matrix, we easily derive from 
and (EToTj) that 



PAP* = E^E;^) ^Sr-x, P 1 U r = U r - 1 . (3.54) 
It is apparent that 

det SV-i ^ 0, PxArP* = A r _i. (3.55) 
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In view of (I3.54p and (I3.55p . the factorization Theorem 4 from (25] (see 
also [27, p. 188]) yields 

w A (r, A) ={i m - ijU^P^PArP* - \i m2 y 1 (ps; 1 p*y 1 ps; 1 n r ^ 

xw A (r-l,X). (3.56) 



Proposition 3.13 The fundamental solution W of the system (ll.ip . where 
W is normalized by the condition (13. ip and the potential {C^} satisfies (II. 2p . 
admits reprezentation 

W r+1 (z) = (l + iz) r+l w A (r,(2z)- 1 ). (3.57) 

Proof. Formulas fll.ip and (12. 9 p imply the following equalities 

W r+1 {z) = (l + iz)(l m + 2iz(l + iz)- 1 n (ry 1 (r))W r (z) (r > 0). (3.58) 

On the other hand, we easily derive from fl3T29|) . fl3T32|) . flgJH]) and f l3~5T|) 
that 

(PA r P* - XL^y 1 = -(A + i/2) _1 J m2 , 57 1 = £ r *£ r , (3.59) 

PSyP* = {e;)*e;, PSyil r = PE*F r = (e")*7(r). (3.60) 

We substitute f l339j) and f |3T60|) into fl336|) to obtain 

iu A (r, A) = (l m + ^A^j 7 ( r )* 7 (r))«; A (r - 1, A) (r > 1). (3.61) 

In a similar way, we rewrite f !3.52p (for the case that r = 0) in the form 

w A (0, A) = I m + ^tJ7(0)*7(0). (3.62) 

Finally, we compare (" 13 . 58[) with (I3.6ip and (I3.62j) (and take into account 
(13. ip ) to see that W\(z) = (1 + iz)wa(0, (2z) _1 ) and iterative relations for 
the left- and right-hand sides of (I3.57P ) coincide. ■ 
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4 Weyl theory: inverse problems 



The values of if and its derivatives at z = i will be of interest in this section. 
Therefore, using (13. 9p we rewrite (13. 2 j) in the form 

p r (2,P) = -[0 / m2 ]iy r+1 (2)*P(^)([/ mi OjWr+a^)^^))" 1 , (4.1) 

where P in (14. ip differs from P in (13. 2 p by the factor j (and so this V is also 
a nonsingular matrix function with property-j). 

Definition 4.1 Weyl functions of Dime system (II. ip (which is given on the 
interval < k < r and satisfies (II. 2p ) are m 2 x m x matrix functions ip(z) 
of the form (14. ip . where V are nonsingular matrix functions with property- j 
such that V{i) are well-defined and nonsingular. 

It is apparent that (14. ip is equivalent to 

= jW r+1 {z)*V(z)([l mi 0]W r+l (z)*V(z)y 1 . (4.2) 



Lemma 4.2 LetV satisfy conditions from Definition ^. 1\ Then we have the 
inequality 

det ( [I mi 0] W P+1 (-i)*P(i)) ± 0. (4.3) 
Proof. First note that in view of (12. 8p we obtain 

J m + C kJ = 2P(k)*P(k)j. (4.4) 
Formulas ( \3A\t and (14 .4ft imply 

[i mi o]w r+1 (-iyv(i)=r +1 ([i mi o]fi(oy) PP)')- 

x (/3(r-l)j/3(r)*)(/3(r)jT(i)). (4.5) 

Using Proposition 13. 101 (and the second equality in ( 12. 8p ) and putting, corre- 
spondingly, i? = /?(&)* and = j5(k + 1)* or $ = 'P(i), we derive inequalities 

det(/3(fc)j/3(fc + 1)*) ^ and det(/3(r)jT(i)) ^ 0, (4.6) 
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respectively. In the same way we obtain det ( [l mi 0] 0(0)*) ^ 0. Now, 
inequality ( 14. 3 p follows from (14. 5p . ■ 

Our next proposition is proved similar to Corollary 13.51 

Proposition 4.3 Suppose if is a Weyl function of Dirac system ( II. ip on 
the interval < k < r, where the potential {C k } satisfies (II. 2p . Then if is a 
Weyl function of the same system on all the intervals < k < r (r < r) . 

Proof. Clearly, it suffices to show that the statement of the proposition 
holds for r = r — 1 (if r > 0). That is, in view of Definition 14. 1[ we should 
prove that V[z) := (I m — izC r j)V(z) has property-j, that V{i) is well-defined 
and that the first inequality in (13. 3p written for V at z = i always holds (i.e., 
V(i) is nonsingular) , if only V has these properties. 
Indeed, since we have 



{I m - izC r j)*j{I m - izC r j) = (1 + \z\ 2 )j + i(z - z)jC r j > (1 + \ Z \ 2 )j, 



the matrix function V has property- j. The non-singularity of V(i) 
C r j)V(i) is apparent from (14 . 4 f) and (14. 6p . ■ 



(4.7) 

(Im + 



Theorem 4.4 Suppose if is a Weyl function of Dirac system ( II. ip on the 
interval < k < r, where the potential {C k } satisfies (II. 2p . Then {C k } k=0 
is uniquely recovered from the first r + 1 Taylor coefficients of <p (ijrf ) oi 
z = 0. 

If ip (ifef) = X]fc=o + 0( zT+1 )> ^en matrices $^,1 are recovered via 
the formula 



$fc,i 



Jo + 01 + • • • + 



(4.8) 



IS 



Using $^,1 we easily recover consecutively H k = [$ fc l $^,2] [where $^,2 i 
given m (I3.34p ) and Sk, which is the unique solution of the matrix identity 
A k Sk — S k A* k = iH k jH* k . Next, we construct 

1 (ky 1 (k) = IliS k r 1 P*(PS k r 1 PT 1 PS k 1 Tl k , P=[q ... J m2 ]. (4.9) 
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Finally, we use r ){k)* r y{k) to recover via ( 12. 9p . 
Proof. Put 

A{z) := |1 + z 2 \- 2{r+1) [I mi <f{z)*\ W r+1 (z)*jW r+1 ( 
According to ( 13. 9 p and (14. 2p we have 



I mi 



(4.10) 



A(z) 



x 



[I mi 0]W r+1 (z)*P(z) 

[i mi o]w r+1 (zyv( 



V(z)*jV(z) 



(4.11) 



From ( 14. 3 p and ( 14. lip we see that A is bounded in the neighbourhood of 
z = i: 

||^(,2)||=0(1) for z-H. (4.12) 

Let us include into considerations the S-node (corresponding to Dirac sys- 
tem), which is constructed in accordance with Remark f3. 121 Substitute ( I3.57P 
into pnUp to obtain 



^)=((i-i2)(i+i^))- r - l [/ mi <p( Z r] 



(4.13) 



x [3 ~ 



n;(^-ii)"V(^-i/)"V) 



where J = I( r +i) m2 - Here we used the important equality 

w A {r,X)*jw A (r ,A) = j - i(\ -X)U.* r (A* r - \I)- l S;\A r -A/)- 1 ^, (4.14) 

which follows from (133U|) and (I332|) (see, e.g., [2HI25]). 

Notice that S r > 0. Hence, formulas (jHTTTjl . (T4~T2"j) and (JUS]) imply that 



^-^ 7 )" n ' 



0(1) for z ->• i. 



(4.15) 



Using the block representation n r = [$ r ;1 $ r 2 ] from ( I3.5ip and multiplying 



both sides of (I4.15P by 



we rewrite the result: 



2z 



2z 



O^ ($; 2 ^ r -^/) _1 $ r , 2 ) 1 ) for 



(4.16) 
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In order to obtain (I4.16P we applied also the matrix (operator) norm inequal- 
ity ||XiX 2 || < ||Xi||||X 2 ||. 

The resolvent (A — XI)~ l is easily constructed explicitly (see, for instance, 
formula (1.10) in [22]). In particular, we derive 

M^-rA^-rk^ ^ ■■■ '»■]■ ?:= tt§'™' 

(4.17) 

From (I4.17P we see that 

K2 ( Ar -±iy\ 2 =i(i-(i^y +1 )i m . (4.i8) 

Partitioning $ rl into m2xm 1 blocks $ r ,i(fc) and using fl4.16p - fl4.18p we obtain 

^f^+^rE^w^n for z ^o, 

^ ' k=0 

which can be easily transformed into 

V (i\^]+(l-z)Y,z k ^ 1 (k) = 0(z r+1 ) for z^O, (4.19) 

and (14. 8 p follows for k = r. Since a(A r ) Ha(A*) = the matrix S r is uniquely 
recovered from the matrix identity (13 .50 p . Finally, (14. 9 j) for the case, where 
= r, is apparent from (I3.60p . From Proposition 14.31 we see that (p is a 
Weyl function of our Dirac system on all the intervals < k < r (r < r) and 
so all CV are recovered in the same way as C r . ■ 

The next corollary is a discrete version of Borg-Marchenko-type uniqueness 
theorems. The active study of such theorems was triggered by the seminal 
papers by F. Gesztesy and B. Simon p^[T9]. 

Corollary 4.5 Suppose (p and <p are Weyl functions of two Dirac systems 
with potentials {Ck} and {Ck}, which are given on the intervals < k < r 
and < k < r , respectively. We suppose that matrices {Ck} and {Ck} are 
positive and j -unitary. Moreover, we assume that 

(p (*Y+ — J ~ ^ y Y+ — / = ^(^ P+1 )' z ~ * °°> P e P — m i n ( r ' 

(4.20) 
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Then we have Ck = Ck for all < k < p. 

Proof. According to Proposition 14.31 both functions <p and (p are Weyl 
functions of the corresponding Dirac systems on the same interval [0, p]. 
From (I4.20p we see that the first p + 1 Taylor coefficients of p> (ifqpj) and 
p> (ifef ) coincide. Hence, the uniqueness of the recovery of the potential 
from Taylor coefficients in Theorem 14.41 yields Ck = Ck (0 < k < p). ■ 

Taking into account (14.81) . we derive that the first r + 1 Taylor coefficients 
of p r (ijzf) at z = (for any Weyl function p r of a fixed Dirac system) can 
be uniquely and in the same way recovered from the matrix $ r i, which, in 
turn, can be constructed as proposed in Remark 13.121 Therefore, the next 
theorem is apparent. 

Theorem 4.6 Let Dirac system (ll.ip . where matrices Ck satisfy (11.21) . be 
given on the interval < k < r . Then all the functions (pd{z) = p r f^i^f > > 
where p> r are Weyl functions of this Dirac system, are non- expansive in the 
unit disk and have the same first r + 1 Taylor coefficients {<pk}o at z = 0. 

Step 1 in the proof of Theorem 13.81 shows that the Weyl function of 
Dirac system on the semi-axis can be constructed as a uniform limit of Weyl 
functions p> r on increasing intervals. Hence, using Theorem 14.61 we obtain 
the following corollary. 

Corollary 4.7 Let (p(z) be the Weyl function of some Dirac system (II. ip . 
which is given on the semi-axis and satisfies (II ,2p . Assume that p r is a Weyl 
function of the same system on the finite interval < k < r. Then the first 
r + 1 Taylor coefficients of <p (i-jrf ) ond <p r (i|rf) coincide. Therefore, the 
system can be uniquely recovered from ip via procedure from Theorem \4-4 

5 Operator identities and interpolation prob- 
lems 

One can easily derive (see, e.g, [T71 p. 474]) that the equality 

Sfc+i,p+i — Sk P = Qkp + Qk+i,p+i — Qk+i,p — Qk,p+i, — 1 < k,p < r — 1 

(5.1) 
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holds for the blocks Sk p and Qk p of the block matrices S r = {sfc p }£ p=0 and 
Qr = {Qkp}k, P =oi respectively, which satisfy the operator identity 

A r S r - S r A* r + \Q = 0, (5.2) 

where A r is given by (I3.32p . Here we add sometimes commas between the 
indices of blocks and put also 

s -i,p = s k,-i — Q-i,p — Qk,-i = 0- (5-3) 

For the case that S r corresponds to Dirac system, we rewrite (15 .11) below in 
an equivalent form and obtain the structure of S r . 

Proposition 5.1 Let S r satisfy (I3.50P . where A r , $ r i and $ r 2 are given by 
(I3.32p . (14. 8 p and the last equality in (I3.34p . respectively. Then S r has the 
following structure: 

SQO = I m2 ~ 0000 aIld S k+l,p+l ~ S k p = 0fc+10* +1 (5.4) 

for — 1 < k,p < r — 1, k + p + 2>0. 

The following statement is immediate from Theorem 14. 61 and Proposition 15.11 

Theorem 5.2 Let Dirac system (II. ip . where matrices satisfy (II. 2p . be 
given on the interval < k < r. Then all the functions tpd{ z ) — fr ^T^f > "Pj ; 
where (p r are given by (13. 2p . matrix functions V(z) in (13. 2 p have property-j 
and matrices V(i) are non-singular, are non- expansive in the unit disk and 
have the same first r + 1 Taylor coefficients {<ftk}o at z = 0. The matrix S r 
determined by these coefficients via (15. 4p is positive. 

On the other hand, if we assume only that the coefficients {4>k}o are fixed 
and S r given (15 .4p is positive, two related interpolation problems appear. 

Interpolation problem I. Describe all the analytic and non-expansive 
in the unit disk matrix functions (fd such that the coefficients {4>k}o are their 
first r + 1 Taylor coefficients. 

Interpolation problem II. Describe all the positive continuations of 
S r , which preserve the structure given by (15 .4p . 
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